Magneto-electrostatic trapping of ground state OH molecules 
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We report the magnetic confinement of neutral, ground state hydroxyl radicals (OH) at a density 
of ~3xl0 3 cm -3 and temperature of ~30 mK. An adjustable electric held of sufficient magnitude to 
polarize the OH is superimposed on the trap in either a quadrupole or homogenous held geometry. 
The OH is confined by an overall potential established via molecular state mixing induced by the 
combined electric and magnetic helds acting on the molecule's electric dipole and magnetic dipole 
moments, respectively. An effective molecular Hamiltonian including Stark and Zeeman terms 
has been constructed to describe single molecule dynamics inside the trap. Monte Carlo simulation 
using this Hamiltonian accurately models the observed trap dynamics in various trap conhgurations. 
Confinement of cold polar molecules in a magnetic trap, leaving large, adjustable electric fields for 
control, is an important step towards the study of low energy dipole-dipole collisions. 

PACS numbers: 33.80.Ps, 33.55.Be, 39.10.+J, 33.15.Kr 



Cold and ultracold molecules are currently subject to 
intense research efforts as they promise to lead major ad- 
vances in precision measurement, quantum control, and 
cold chemistry [lj. Much of this interest has focused 
on molecules possessing permanent electric dipole mo- 
ments. The long range, anisotropic dipole-dipole interac- 
tion leads to novel collision 0, 0] and reaction [1] dynam- 
ics, and may serve as the interaction between molecular 
qubits [5|L in architectures for quantum information pro- 
cessing [6| . Quantum phase transitions can be explored 
with polar molecules [7|. Furthermore, polar molecules 
may possess extremely large internal electric fields (on 
the order of GV/cm), yielding dramatic sensitivity en- 
hancement to searches for an electron electric dipole mo- 
ment [1, Q . Cold polar molecules have been produced 



via photoassociation of ultracold atomic species llOlllll , 
Stark deceleration , [13] , and buffer gas coolingjl^. 

Polar molecules are readily trapped in inhomogeneous 
electric fields [ljl however theory suggests that 
at high phase-space densities such traps invariably suf- 
fer from large inelastic collision losses [Tjj. One pos- 
sible solution is to trap strong-field seeking states via 
time-dependent electric traps |l8|. In contrast, we are 
interested in magnetically trapping these neutral polar 
molecules so as to have the freedom to apply external 
electric fields to control collision dynamics inside the 
trap. Highly vibrationally excited, triplet KRb molecules 
produced via photoassociation have been magnetically 
trapped at a temperature of ^300 /j,K and density of 
10 4 cm" 3 [19]. Buffer gas cooling has been used to 
load a magnetic trap with ground-state CaH and NH 
molecules at temperatures of hundreds of mK and den- 
sities of ~10 8 cm -3 , but in the presence of a cold He 
buffer gas [13, 0] . We report here the first observation 
of magnetically trapped hydroxyl radical (OH) molecules 
and the lowest temperature (30 mK) yet achieved for 
a magnetically-trapped polar molecule in its rovibronic 
ground state. Moreover, we perform the first study of 




FIG. 1: (color online) Illustration of the magneto-electrostatic 
trap (MET) assembly. The terminus of the Stark decel- 
erator (yellow, gray electrodes), shown to the left, couples 
state-selected cold molecules into the magnetic quadrupole 
(blue). The double-electrode structures within the electric 
quadrupole (red) allow for application of uniform electric 
fields within the magnetic trap. The ring (green) supports 
a 1 inch lens for collection of laser-induced fluorescence. 



trap dynamics based on a single molecule possessing large 
magnetic and electric dipole moments in B- and E-fields 
that are inhomogeneous and anisotropic. Interesting po- 
lar molecule dynamics have been predicted at large field 
strengths [2l[ . Accurate correspondence between experi- 
mental data and simulations is achieved only by account- 
ing for the molecular state mixing induced by the crossed 
fields. Large magnetic (B) fields may serve to suppress 
inelastic collision losses [22J while at the same time sym- 
pathetic cooling of molecules via co-trapped ultracold 
atoms may become feasible 23[. In addition, the mag- 
netic quadrupole trap described here permits the appli- 
cation of arbitrary external electric (E) fields to a large 
class of polar molecules — a necessary step towards ob- 
serving cold molecule dipole-dipole collisions 24 1 . 
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To magnetically trap OH, we begin with a cold beam 
of ground state ( 2 n 3 / 2 ) OH produced via a pulsed elec- 
tric discharge through 2% H 2 vapor [25J] seeded in 1.5 
bar of Kr. A piezoelectric transducer valve [2(| operating 
at 10 Hz provides a supersonic expansion of the OH/Kr 
mixture through a 1 mm nozzle. This results in a 490 
m/s molecular packet with a 15% longitudinal (z) ve- 
locity spread. The transverse (p) velocity spread of the 
beam is limited to 2% by the 3 mm diameter molecular 
skimmer. Compared with a Xe expansion, the higher- 
velocity Kr expansion leads to less transverse beam loss 
within the current decelerator design. Thus, despite the 
fact that the Kr beam requires a higher phase angle for 
deceleration (smaller decelerator acceptance) it results in 
more decelerated molecules. After passing through the 
skimmer, the OH packet is spatially matched to the de- 
celerator entrance via an electrostatic hexapole. Details 
regarding the OH Stark decelerator design and operating 
principle may be found in previous work 27 , 28| . We uti- 
lize the 142 stages of our Stark decelerator to slow a 130 
mK packet (in the co-moving frame) of OH to 20 m/s at 
a phase angle <j>o = 47.45° for coupling into the magneto- 
electrostatic trap (MET). Decelerated OH packets are 
detected via laser-induced fluorescence (LIF) upon en- 
tering the MET, whose center lies 2.7 cm from the final 
pair of decelerator rods. The 282 nm pulsed laser beam 
is oriented along the mutual longitudinal axis (z) of the 
decelerator and MET, and the resulting 313 nm LIF is 
collected in a solid angle of 0.016 sr. 

The MET is illustrated in Fig. [IJ and consists of two 
copper coils operated in an anti-Hclmholtz configuration 
surrounded by an electric quadrupole. The center-to- 
center magnetic coil spacing is 1.5 cm, providing a field 
gradient of 6700 G/cm by passing 1500 A through the 4 
turns of each water-cooled coil. A longitudinal trap depth 
of 400 mK is achieved, while the electrodes surround- 
ing the coils allow the application of an E-field to the 
trapped sample. The double- and single-rod electrodes 
that comprise the electric quadrupole allow the applica- 
tion of two distinct field configurations. A uniform bias 
electric field is created by grounding the single-rods while 
oppositely charging the opposing double-rods to ±3 kV. 
Alternatively, we generate a quadrupole E-field by charg- 
ing neighboring rods with opposite polarities of ±9 kV. 
Shielding by the MET coils reduces the interior E-field 
by a factor of ~5. Figure [2] depicts the quadrupole elec- 
tric field superimposed on the magnetic quadrupole field. 
OH has both an appreciable magnetic and electric dipole, 
and we can easily alter the geometry of the trapping po- 
tential by applying electric fields of variable magnitude 
and direction to the magnetically trapped sample. 

The OH ground state is best described by Hund's case 
(a), in which the spin (|S| = i) and orbital (|A| = 1) 
angular momenta are strongly coupled to the molecu- 
lar axis, yielding the total angular momentum projection 
Ifil = |. The Ifil = i state lies 126 cm" 1 above this 



a), 



& 0.4 
f 0.2 

£ o.o 

!-»-* 

-0.4 

-o.e 

-0.B 







0.5 




^^^^-m = 312 


0.4 
0.3 






g 0.2 






i o.i 






m 






■C Q.O 




a 






<fl 






-0,2 






-0.3 






-0,4 





FIG. 2: (color online) (a) Zeeman and (b) Stark effects of the 
ground state structure of OH. (c) Magnetic and (d) electric 
quadruple fields viewed from z. (e) Side view of the MET 
configuration with quadrupole E and B fields, (f) Transverse 
adiabatic potential surfaces for various components of the OH 
ground state at the longitudinal MET center. The top surface 
depicts the decelerated/trapped |J = §,mj = §) state. 



ground state. The nuclear spin of the hydrogen atom 
(I = |) leads to a hyperfine splitting of the ground 
state into F = 2, 1 components. In the presence of large 
external fields (E>1 kV/cm or B>100 G), the electron 
angular momentum and nuclear spin decouple and 
one considers the total angular momentum J, which in- 
cludes nuclear and electronic orbital angular momentum 
as well as electron spin. The projections mj of J are 
defined along the axis of the applied field. The ground 
state is further split into two opposite-parity states (//e) 
by the coupling of electronic orbital angular momentum 
to nuclear rotation. This small A-doublet splitting of 
1.67 GHz, along with the 1.67 D electric dipole, is re- 
sponsible for the rather large Stark shift experienced by 
OH molecules [29| . Fig. HJa) and (b) show the Zeeman 
and Stark effects of the OH ground-state respectively. 
External electric fields can thus dramatically alter the 
potential confining the magnetically trapped molecules. 
Fig.[Hc-e) illustrate the field geometries within the MET. 
In the plane perpendicular to z, Fig. [^c) shows the ra- 
dial B-field characteristic of a magnetic quadrupole, while 
Fig. IHd) depicts the electric quadrupole field introduced 
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by the electrode geometry. Fig. [2^e) illustrates both B 
(blue) and E (red) fields along z. 

The top potential surface of Fig. [5£f) represents the 
energy shift of the decelerated \J = §, raj = |) state (E- 
fields couple e and /, therefore they are no longer good 
quantum numbers) within the combined quadrupole E- 
and B-fields present in the MET. The transverse profiles 
of the four trapping adiabatic potentials are depicted at 
the longitudinal trap center, where B ■ z = 0. The re- 
maining four surfaces (not shown) are simply inverted rel- 
ative to those displayed, and are therefore anti-trapping. 
The three-dimensional MET potentials are calculated by 
diagonalizing an effective Hamiltonian that includes both 
Stark and Zeeman terms. Coupling between the ground 
| J = |,0 = §}, and excited |§, §), ||,|), ||, |) states 
is included, for a total of 64 hyperfine components. The 
| J = |,fi = |), |§, 3) coupling is the most critical, in- 
creasing the trap depth of the top surface of Fig. 0(f) 
by 11%. Addition of the remaining states modifies the 
potential by <0.1%. The square shape of the top surface 
is a result of the varying angle between E and B in the 
MET. The potentials from each field directly sum where 
the fields are collinear, but the trap shallows from this 
maximum value as the angle between E and B increases. 

The Stark decelerator slows only those OH molecules 
in the electrically weak-field seeking (EWFS) component 
of the ground state, i.e., the upper A-doublet. Of these 
EWFS molecules, only those with components J = |, 
rrij = ±| are synchronously slowed to 20 m/s. Because 
the positive mj state is magnetically weak-field seeking, 
only 50% of the decelerated molecular packet is trappable 
by the MET. The molecules experience a B-ficld >100 G 
at the last deceleration stage emanating from the back 
coil, preserving the quantization of the decelerated state. 

To trap the incoming OH packet, the coil closer to the 
decelerator (front coil) is grounded while the coil further 
from the decelerator (back coil) operates at 2000 A for 
the duration of the 3.7 ms deceleration sequence. The 
effective magnetic dipole of OH in its J = mj — | state 
is 1.2 hb — 0.81 K/T, where /i^ is the Bohr magneton, 
allowing the back coil to stop molecules with velocity 
< 22 m/s. The front coil is turned on after the 20 m/s 
molecules are stopped by the back coil, which occurs 2.65 
ms after exiting the final deceleration stage. The trap 
switching is synchronous to the 20 m/s molecules since 
OH possessing larger longitudinal velocities escape the 
quadrupole field. For technical reasons, the two magnet 
coils are connected in series, causing the current in both 
coils to decrease from 2000 A to 1500 A over ~800 fxs once 
the front coil is switched into the circuit. This problem 
will be addressed with a current servo. 

Figure [3] displays typical time-of-flight data result- 
ing from the trapping sequence along with results from 
Monte Carlo simulation of the trap dynamics. The large 
initial peak is the decelerated OH packet traversing the 
detection region at the center of the trap. All subsequent 
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FIG. 3: (color online) Time-of-flight data and Monte Carlo 
simulation results (solid red line) for two different MET con- 
figurations. The magnetic field switches from the stopping 
configuration to quadrupole trapping at t = 2.65 ms. (a) Mag- 
netic trap only, (b) Stopping and trapping in the presence of 
combined electric and magnetic quadrupole fields. Note larger 
steady-state trap population. 



peaks are the result of oscillation in the p dimension, with 
a frequency of ~650 Hz. Fluorescence is gathered in p, 
and oscillations in z are not detectable since our detection 
region spans the entire visible trap length. The magnetic 
quadrupole trap is unable to confine all the molecules 
stopped by the back coil alone due to the 25% current 
drop, as is clear from the difference in signal height at 
~2.5 ms versus the steady-state level at, e.g., 10 ms. 
Note the larger steady-state population of Fig. [3^b) from 
the addition of a confining transverse electric quadrupole. 
Also visible is the shorter coherence time compared to 
Fig. [UJa), which results from trap distortion visible in 
Fig. [H(f). An estimate of the density of trapped OH, as 
measured by LIF, gives ~ 3xl0 3 cm -3 . By comparison, 
the density of a packet "bunched" at a phase angle of 
0o = 0° is 10 7 cm" 3 , while the 20 m/s packet is 10 4 
cm" 3 . Monte Carlo analysis of the velocity distribution 
of the trapped OH yields temperatures of ~30 mK. With 
our improved knowledge of the trap dynamics due to ac- 
curate Monte Carlo simulations, we expect to be able to 
trap 2-3 times more OH through modification of decel- 
eration phase angle and trap switching times for optimal 
coupling of the slowed beam into the MET. 

Higher decelerator phase angles should, in principle, 
produce a slower OH packet. However, below 20 m/s the 
OH population drops sharply and no further deceleration 
is observed. We attribute this effect to a reduction in the 
effective aperture seen by slow molecules as they traverse 
the final deceleration stage. The OH closest to the final 
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the different surfaces depicted in Fig. [2jf) akin to Majo- 
rana transitions observed in atomic magnetic quadrupole 
traps. Furthermore, higher densities within the MET will 
facilitate the study of cold OH-OH collisions in variable 
electric fields. As rather modest electric fields are ex- 
pected to modify collision cross sections between polar 
molecules by as much as 10 3 17|, crossed-beam exper- 
iments exploiting the MET described here will benefit 
from both low center-of-mass collision energies and tun- 
able E-fields in the interaction region. 

We acknowledge DOE, NIST, and NSF for support. 
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technical assistance. B. Lev is a NRC postdoctoral fellow. 



FIG. 4: (color online) Measured (r) and de-convolved colli- 
sional (r c ) lifetimes of magnetically trapped OH at a back- 
ground pressure of 1 x 1CP 6 Torr (•, red line) and 4 x 1CP 8 
Torr (o, blue line). 



decelerator rods see a larger longitudinal potential than 
those on-axis, and at sufficient phase angles, are stopped 
or reflected and hence cannot be observed in the MET re- 
gion. Numerical models elucidate this effect, which may 
be mitigated by new decelerator designs [301 ] . 

The 10 Hz repetition rate of our Stark decelerator cur- 
rently limits our trap interrogation time to 100 ms. Main- 
taining continuous currents of 1500 A presents technical 
challenges due to heating in both the switching tran- 
sistors and cabling. To mitigate this, we operate the 
switching MOSFETs on a water-cooled copper plate and 
also run chilled water through the MET coils. Water- 
cooled cables are also necessary for such high-power oper- 
ation. Figure 2] shows electric field- free OH trap lifetimes 
measured in two different background pressure regimes: 
1 x 10~ 6 and 4 x 10~ 8 Torr, both dominated by N 2 at 
295 K. Non-ideal power supply control causes the trap 
current to decrease from 1500 A to 1160 A at a rate of 
2000 A/s. This lowers the trap depth, resulting in a loss 
of 12 ± 2 s _1 . After accounting for this loss mechanism, 
we find background collision-limited lifetimes of 20 ± 5 
ms and 500 ± 130 ms for the two pressure regimes, re- 
spectively. These lifetimes determine an OH-N 2 collision 
cross section of 500 ± 100 A 2 . From this clear dependence 
of trap lifetime on background gas pressure, we envision 
measuring dipole-dipole collision cross sections by using 
polar molecules as background and polarizing them with 
the E-field of the MET. 

We have magnetically confined ground state hydroxyl 
radicals in a MET that allows for the application of a 
variable E-field to the trapped polar molecules. Diago- 
nalization of an effective Hamiltonian modelling the en- 
ergy shift of OH in combined E- and B-fields was nec- 
essary for accurate Monte Carlo simulation of the trap 
dynamics. Future experiments could quantify temper- 
ature dependent non-adiabatic transition rates between 
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